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We present a general formulation of the nonlinear, time-dependent current 
through a small interacting region, where electron energies are changed by time- 
dependent voltages. An exact solution is obtained for the non-interacting case when 
the elastic coupling to the leads is independent of energy. Temporal phase coherence 
in a double-barrier tunneling structure produces "ringing" in the response of the cur- 
rent to a voltage pulse, which can be observed experimentally in the dc-current by 
varying the pulse length in a train of voltage pulses. The nonlinear current due to 
an ac-bias also shows complex time-dependence. 

PACS numbers: 73.20.Dx 73.40.Ei 73.40. Gk 73.50.Fq 



Typeset Using REVTEX 



1 



The importance of the spatial coherence of electronic wavefunctions is one of the hall- 
marks of the mesoscopic regime. A panoply of mesoscale phenomena including weak- 
localization and weak anti-localization, 1 Aharonov-Bohm oscillations, 2 and Universal con- 
ductance fluctuations 2 all rely on the phase coherence of electrons in small structures. Since 
the traditional experimental probe of these effects is steady-state transport, the role of tempo- 
ral phase coherence is generally subsumed under the effects of spatial coherence. However, 
recent experimental progress 3-5 has opened the door to direct measurement of electronic 
phase coherence in time. 

Interest in the time domain has sprung both from device applications of double-barrier 
resonant-tunneling structures, 3 and from possible current standard applications of single- 
electron tunneling circuits. 4 While the interactions among electrons are mean-field like in 
the three-dimensional resonant-tunneling structures, the low-dimensional confinement of 
electrons in the single-electron circuits results in strong many-body correlations. A general 
theory for temporal coherence in these structures must therefore rigorously account for 
interactions. 

In this letter, we first present a general formulation of the nonlinear, time- dependent 
current through a small interacting region coupled to two non-interacting leads. 6 A general 
result, Eq. (g), expresses the time-dependent current flowing into the interacting region from 
one lead in terms of local Green functions. This expression is then used to explore temporal 
coherence in the response of a mesoscopic system to time-dependent external driving. We 
focus on a fully nonlinear, exactly solvable example — non-interacting electrons traversing 
a double-barrier tunneling structure with energy independent coupling to the leads. The 
temporal coherence of this system is evident in the "ringing" of the current in response to a 
rectangular pulse of the bias. This ringing can be observed experimentally in the dc-current 
by varying the pulse length in a train of voltage pulses applied to the structure. Similarly, 
the current flowing in response to an ac-bias displays complex time-dependence, which is 
reflected in oscillations of the dc-current vs. driving frequency. 7 

The formalism and Hamiltonian are similar to those employed for the steady-state current 
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in Ref. 8, 

H= £ka(tH a C ka + H int [{d n },{dl},t] + [VkaAtHadn + h.C.], (1) 

k,a&L/R k,a£L/R 

n 

where c^ a {c ka ) creates (destroys) an electron with momentum k in channel a in either the 
left (L) or the right (R) lead, and {d+} and {d n } form a complete, orthonormal set of 
single-electron creation and annihilation operators in the interacting region. However, time 
dependence due to external driving is now explicitly included in the energies of states in the 
interacting region and in the leads, in the hopping matrix elements, and in the interactions 
themselves. A central assumption is that interactions between electrons in the leads and 
between electrons in the leads and in the interacting region can be neglected. Geometrically, 
the leads must therefore rapidly broaden into large metallic contacts in which interactions 
are strongly screened. Since transport experiments on mesoscale structures typically satisfy 
this requirement, the model is directly relevant to experiment. Note that because time- 
dependent potentials must be associated with changes of charge, there will be capacitive 
currents flowing in the contacts in addition to the currents from the interacting region 
determined by ([!]). The contribution of the capacitive currents can be eliminated, however, 
by measuring only the time- averaged current. 

The Keldysh approach 9 to calculating the current flowing into and out of the interact- 
ing region treats the contacts as systems separately in equilibrium in the distant past, 
possibly with different chemical potentials. Physically, applying a time-dependent bias 
(electrostatic-potential difference) changes the energies of states in the leads via efc Q (£), with- 
out changing their occupations. This preserves the coherent evolution of phase in the leads, 
^ka(t) oc exp[— i J dt' ek a (t')}, and produces interference in tunneling between the leads and 
the interacting region, due to their different time- dependent energies. 

The time-dependent current from the left lead into the interacting region is 

J L (t) = ^Re{ E V ka , n (t)G< ka (t,t)}, (2) 

71 

where in (§) we have just rewritten the expectation of the current matrix element in terms of 
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the Keldysh Green function 9 G^ ka (t,t') = i(c ka (t')d n (t)) . Since the Hamiltonian describing 
the leads is non-interacting, one has the Dyson equation 

^n,fca(*) *') — 2J / ^l^to,m(*l)[^,m(Ml)#fca,ta^ (3) 
m 

where G< m (t,t') = z(d+ ,(f)d n (t)), and G r n>m (t,t') = -i6(t - i?){{d+ m {i?),d n {t)}) is the re- 
tarded Green function. The time-dependent Green functions in the leads for the uncoupled 
system, which appear in (Q), are given by 

t 

9L,ka^ tn >^ i6 ^' '-t)exp\-if dtie fcQ (ti)], (4) 

L Jt' J 

where is the energy of the state k, with occupation /(e° a ), when the system was prepared 
in the distant past, and the advanced Green function is related to the retarded Green function 
by g a (t,t') = [g T (t',t)\*. It is convenient in (0) to turn the sum over momentum states k in 
the leads into an integral over energies, and to define the elastic coupling between the leads 
and the states in the interacting region via 

[r L (e, t', t)] m>n = 2tt £ Pa (e)V a>n (e, t)V£ m (e, t') exp [i f dt 1 A«(e, t x )\ , (5) 

where p a (e) is the density of states in channel a and the energy, e/- a (t), of each state in the 
leads is separated into a constant part e° ka = e and a time-dependent part A Q (e, t). Rewriting 
the current in terms of the elastic coupling, and using matrix notation for the level indices 
in the interacting region, we find 

J L (t) = ~ f dt> [ ^ImTr{e^- t ')r i (e,t',t)[G<(t,t') +/L(e)G r (t,t')]}- (6) 

An analogous expression applies for the current flowing in from the right lead. 

Eq. (Q) is the central formal result of this work. The strong resemblance to the steady- 
state result in Ref. (8) means that the time-dependent problem is not significantly harder 
than the time-independent one. For a range of mesoscopic systems, a rigorous formula 
for the steady-state current is a useful tool, and we believe (0) will be a useful tool to 



study dynamics. While interacting quantum transport is addressable with (H), we use it 
here to explore dynamics in an exactly solvable, non-interacting system corresponding to a 
quantum-well structure. 

In general, if interactions are neglected, the retarded and advanced Green functions are 
given by standard Dyson equations, 10 and the Keldysh Green function in (|5]) is related to 
them via 



G < (t,t')=ijdt 1 Jdt 2 G r {tM)[L J ^e^-^f L/R {t)T L l R {eMM)]G a {t2,t'). (7) 



de ic 

A tractable, non-interacting example of particular interest corresponds to a double-barrier 
tunneling structure containing a single resonant level, eo(t) = 6q + A(t), with energy- 
independent coupling to two leads, 



r L/R (t',t) = T L/R expfz/^Ai/^!)]. (8) 



In (H), the time dependence in the leads is restricted to a rigid shift of all states in the 
left (right) lead by A^^it). Within this framework, application of a voltage bias corresponds 
to a shift of energies in one lead with respect to the other, and generally a shift of the resonant 
level as well. For energy-independent coupling to the leads, the retarded Green function for 
the resonant level is independent of energy shifts Ai/R{t) in the leads, and is given by 11 
G r (t,t') = exp[— T(t — t')/2] g r (t,t'), where g r (t,t') is the retarded Green function for the 
uncoupled level, 



g r (t,t') = -i9(t-t') exp[-z /'dtieoOi)!, (9) 

L Jt' J 



and where r = Tl + T R is the total elastic coupling to the leads. 

The resulting expression for the retarded Green function of the resonant level, G r (t,t'), 
can be used in (0) to generate the Keldysh Green function, G < (t,t'), which contains all 
the information concerning non-equilibrium occupations. Physical quantities such as the 
occupancy of the level and the currents through the barriers can then be expressed directly 
in terms of G < (t,t'), or, more conveniently via (0), in terms of G r (t,t') and the occupation 



functions in the leads. The occupancy of the level, n(t) = —iG < (t,t), and the current, Jl/r 
from (|jD, are given by: 

n(t) = £ T L/R [ ^f L/R (e)\A L/R (e,t)f (10) 

L/R J Z7F 

Jl/rW = ~V L/R [n{t) + J ^f L/R {e)lm{A L/R {e,t)}], (11) 
where we have defined 

A L/R (e,t) = f dhexpli f dt 2 (e + A L/jR (t 2 ))lG r (t,ti). (12) 

In the time-independent case, A L / R (e,t) is the Fourier transform of the retarded Green 
function, and one finds the usual result for the resonant-tunneling current, 

J = T fde [h(e) - f R (e)]- 2 - (13) 



hJ ^- J(e _ eo)2 + (r/2)2 

While arbitrary time dependence of both the bias and level energy can be addressed via 
Eqs. ([^-[12|), we have chosen to consider the responses to a rectangular pulse and to an ac- 
bias as examples of experimental relevance. As discussed above, an important experimental 
caveat is that the time-dependent current flowing in the contacts will include capacitive 
contributions from the accumulation and depletion layers on either side of the tunneling 
barriers. 12 It is possible, however, to measure exclusively the current flowing through the 
barriers by measuring the time-averaged current, since the dc-current must be uniform 
throughout the structure. 

From Eq. flT2]) , one finds for a rectangular pulse of duration s starting at t = 0, 

e - e + iT/2 - (A - A L/R ) [l - e i(«o-HT/2)(t-«) (l - e i(e-eo-A+A i/H+ ir/2) S 



ALIR ^ t] (e - e + iV/2) (e - e - A + A L/R + iV/2) 

(14) 

Expression ( |i~4"D for A L / R (e, t) applies for times, t > s, after the pulse has ended. For times 
during the pulse, < t < s, one should replace s by t, which follows from causality. In 
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Fig. la, the current flowing through the barriers is plotted for a voltage pulse of duration 
s = 3h/T (dashed curve). Before the pulse, the chemical potentials Hl/r an d the level energy 
eo are equal so the current is zero. During the pulse, energies in the left lead are increased 
by Al = 10r and the energy of the resonant level is raised by A = 5r, appropriate for a 
symmetric structure (inset, Fig. lb). (The choice A R = is arbitrary since only the relative 
shifts of energy are significant.) Initially the current through each barrier grows proportional 
to £log(l/t), 13 and then oscillates with a period At = 2irh/\(^ L / R + A L / R ) — (e + A)| (~ 1.26 
in Fig. la). The time scale for the decay of the oscillations is the resonance lifetime h/T, 
i. e., the time during which the states comprising the resonance remain in phase and hence 
contribute constructively to the current oscillations. For the case shown in Fig. la, the 
occupancy remains fixed at n = 0.5 by symmetry, and so the currents through the two 
barriers must be equal. The "ringing" in the current is a consequence of the different phase 
evolution of the resonant level and of the states in the leads. Such phase coherence is 
explicitly absent in the calculation of time-dependent resonant-tunneling currents by Chen 
and Ting, 14 and so they find only exponential relaxation following an abrupt change of bias. 

To observe the "ringing" experimentally, a series of pulses such as that of Fig. la could 
be applied to a tunneling structure and the dc-current measured as a function of pulse 
duration. 15 In Fig. lb, the derivative of the dc-current with respect to pulse length is 
plotted, normalized by the repeat time, r, between pulses. For pulse lengths, s, of the order 
of the resonance lifetime, h/T, the derivative of the dc-current mimics closely the time- 
dependent current following the pulse, and, likewise, asymptotes to the steady-state current 
at the new voltage. We assume that the pulse duration, of order a typical resonance lifetime, 
h/T ~ 1 — lOOOps, will be much shorter than the time between pulses, so that pulses will 
not interfere. 

Instead of applying a bias pulse, it may be experimentally more practical to apply an 
ac-bias to the tunneling structure. 3 From Eq. (0), one finds for an ac-potential, 

A L/R {e,t) = exp h—f-^ * (— e _ eo _ knu + ,t /2 » ( 15 ) 
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where A L / R {t) = A L / R cos(ut), A(t) = Acos(ut), and Jk is the k th order Bessel function. 
In Fig. 2a, the current is plotted for an ac-potential of frequency uj = 2T/h. While the 
current has the same period, 2tc/uj, the complex time- dependence inside each period is 
similar to the "ringing" response to a pulse shown in Fig. 1. One consequence of this 
complex harmonic structure is that for temperatures k B T < hu the time-averaged current, 
Jdc, oscillates as a function of period, as shown in Fig. 2b. An oscillation occurs whenever a 
photon-assisted resonant-tunneling peak aligns with one chemical potential, eo + khu = Hl/r, 
giving oscillations periodic in 1/u. 7 Direct measurement of the ac-current is complicated by 
the capacitive contributions from the leads, but we expect the structure shown in Fig. 2a 
to survive. 

In conclusion, a general formula (]6|) for the time-dependent current through an interact- 
ing mesoscopic region has been obtained in terms of local Green functions. While for the 
interacting case only approximate Green function solutions may be available, for a single 
non-interacting level with energy-independent coupling to two leads an exact solution has 
been obtained - |T2|). This case corresponds experimentally to resonant tunneling through 
a double-barrier structure. We find that temporal coherence of electrons tunneling through 
the resonant level leads to "ringing" of the current in response to an abrupt change of bias, 
and similarly complex time-dependence in response to an ac-bias. This "ringing" can be 
observed experimentally in the dc-current by varying the pulse length in a train of voltage 
pulses. We hope that time-dependence will provide a new window on coherent quantum 
transport and will lead to significant new insights in the future. 
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Figure Captions: 

(la) Time-dependent current, J(t), through a symmetric double-barrier tunneling structure 
in response to a rectangular bias pulse. All energies are in units of the elastic coupling to 
the leads, T, the current is in units of eF/h, and all times are in units of Ti/T. Initially, the 
chemical potentials hl and /ir and the resonant-level energy e are all zero. At t — 0, a bias 
pulse (dashed curve) suddenly increases energies in the left lead by A^ = 10 and increases 
the resonant-level energy by A = 5 (see inset). At t — 3, before the current has settled 
to a new steady value, the pulse ends and the current decays back to zero. By symmetry, 
the occupancy of the resonant level is always 0.5 and the currents through the two barriers 
are equal. The temperature in all figures is k B T = 0.1T. (lb) Derivative of the integrated 
dc-current, Jd c , with respect to pulse duration, s, normalized by the interval between pulses, 
r. For pulse durations much longer than the resonance lifetime h/T, the derivative is just 
the steady-state current at the bias voltage, but for shorter pulses the "ringing" response of 
the current is evident. 

(2a) Time-dependent current, J(t), through a symmetric double-barrier tunneling struc- 
ture for an ac-bias of frequency cu = 2T/h (dashed curve). The ac-driving amplitude is 
A L = 10 about fi L = 10 in the left lead, A = 5 about e = 5 for the level, and A R = fi R = 
in the right lead (see inset). By symmetry n = 0.5 and the currents through the two barriers 
are equal. (2b) Time-averaged current, J^ c , as a function of the ac-oscillation period 2n juj. 
The ac-amplitudes are the same as those used in (a). 
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